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Abstract. We study the spatial distributions of /3 + -activity produced by 
therapeutic beams of 3 He and 12 C ions in various tissue-like materials. The 
calculations were performed within a Monte Carlo model for Heavy-Ion Therapy 
(MCHIT) based on the GEANT4 toolkit. The contributions from positron- 
emitting nuclei with T 1/2 > 10 s, namely lc > n C, 13 N, 14 - ls O, 17 . 18 F and 30 P, 
were calculated and compared with experimental data obtained during and after 
irradiation, where available. Positron emitting nuclei are created by 12 C beam in 
fragmentation reactions of projectile and target nuclei. This leads to a /3+-activity 
profile characterised by a noticeable peak located close to the Bragg peak in the 
corresponding depth-dose distribution. This can be used for dose monitoring in 
carbon-ion therapy of cancer. On the contrary, as the most of positron-emitting 
nuclei are produced by 3 He beam in target fragmentation reactions, the calculated 
total /3 + -activity during or soon after the irradiation period is evenly distributed 
within the projectile range. However, we predict also the presence of 13 N, 14 0, 
I7,i8p crea t e( j m charge-transfer reactions by low-energy 3 He ions close to the 
end of their range in several tissue-like media. The time evolution of /3+-activity 
profiles was investigated for both kinds of beams. We found that due to the 
production of 1S F nuclide the /?+-activity profile measured 2 or 3 hours after 
irradiation with 3 He ions will have a distinct peak correlated with the maximum 
of depth-dose distribution. We also found certain advantages of low-energy 3 He 
beams over low-energy proton beams for reliable PET monitoring during particle 
therapy of shallow located tumours. In this case the distal edge of /3 + -activity 
distribution from 17 F nuclei clearly marks the range of 3 He in tissues. 
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1. Introduction 

Beams of protons and carbon ions are used in particle therapy of deep-seated tumours 
for conformal irradiation of a tumour volume while sparing surrounding healthy tissues 
and organs at risk (Castro et al 2004, Amaldi and Kraft 2005). New facilities for 
proton and ion therapy of cancer are planned or under construction in France (Bajard 
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et al 2004), in Italy (Amaldi 2004), in Austria (Griesmayer and Auberger 2004) and 
in Germany (Haberer et al 2004). Along with other facilities the Heidelberg Ion 
Therapy Center (HIT) (Haberer et al 2004, Heeg et al 2004) will use several beams 
for treatment: protons, 3 He, carbon and oxygen nuclei. There oncologists will have 
a variety of treatment tools at their disposal, while each treatment option will be 
characterised by its specific effectiveness, possible side effects and treatment costs. 
Appropriate quality assurance methods should be also developed specifically for each 
kind of treatment. 

In particular, the Positron Emission Tomography (PET) during or after 
irradiation provides the possibility to monitor the delivered dose. The PET monitoring 
methods in proton and ion therapy can be divided into two categories: (i) based on 
tracing the positron emitting nuclei, e.g. 10 C, n C and 15 0, created by proton beams 
in tissues due to fragmentation of target nuclei] (ii) based on tracing the positron 
emitting nuclei, 10 C and n C, created in fragmentation reactions of 12 C beam nuclei. 
The spatial distribution of positron emitting nuclei is measured by detecting gamma 
pairs from e + e~ — > 77 annihilation events. By comparing the measured /3 + -activity 
distribution with the distribution calculated for the planned dose, one can control the 
accuracy of the actual treatment. 

Following extensive theoretical and experimental studies with carbon beams 
(Enghardt et al 1992, Pawelke et al 1996, Pawelke et al 1997, Ponisch et al 2004, 
Parodi 2004), in-beam PET monitoring is successfully used in carbon- ion therapy at 
GSI, Darmstadt, Germany (Enghardt et al 2004, Schulz-Ertner et al 2004). Similar 
approaches can be used for monitoring of proton therapy, as shown early by Bennett 
et al 1975, Bennett et al 1978 and later by Oelfke et al 1996, Parodi and Enghardt 
2000, Parodi et al 2002, Parodi 2004, Nishio et al 2005, Parodi et al 2007a and Parodi 
et al 2007b. PET images from proton and carbon- ion therapy were also studied in 
experiments by Hishikawa et al 2004 at Hyogo ion therapy centre in Japan (Hishikawa 
et al 2002). 

Beams of nuclei lighter than carbon, e.g. 3 He, 4 He or 7 Li, are also of clinical 
interest. This is shown, in particular, by Furusawa et al 2000 and Kempe et al 2007. 
An advantage of 3 He nuclei consists in their specific Z/A = 2/3 ratio which helps to 
protect the 3 He beam from contamination with 4 He, 12 C, 16 nuclei. The feasibility 
of in-beam PET for 3 He therapy was demonstrated for the first time in experiments 
by Fiedler et al 2006. 

In the present work we use a Monte Carlo model for Heavy-Ion Therapy (MCHIT) 
(Pshenichnov et al 2005, 2006) based on the GEANT4 simulation toolkit (Agostinelli 
et al 2003, Allison et al 2006) to study the /3 + -activity profiles induced by 3 He and 12 C 
beams in tissue-like media. We argue that specific nuclear reactions, namely proton 
pick-up by target nuclei, play certain role in production of positron emitting nuclei 
by 3 He beams in addition to previously studied nuclear fragmentation reactions. In 
Section [5] we describe the physical models from the GEANT4 toolkit used to build 
MCHIT. In Section [3] the time-dependent analysis of the /3 + -activity distributions 
induced by 3 He and 12 C beams in graphite, water and PMMA phantoms is presented. 
In Section |4] calculational results are compared with experimental data obtained by 
Fiedler et al 2006. Results for homogeneous phantoms with stoichiometric composition 
of muscle and bone tissues are presented in Section The calculated distributions 
of /3 + -activity induced by low energy proton, 3 He and 12 C beams are discussed in 
Section [6] with emphasis on the role of proton pick-up reactions induced by 3 He. 
In Section [7] the reliability of MCHIT results is verified by comparison with available 
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experimental data on specific reaction cross sections and isotope yields in thick targets. 
Section [5] contains summary and conclusions. 

2. GEANT4 physics models used in MCHIT 

We have used the version 8.2 of the GEANT4 toolkit (GEANT4-Webpage 2006) to 
build a Monte Carlo model for Heavy-Ion Therapy (MCHIT). The model is intended 
for calculating the spatial distributions of dose and /3 + -activity from beams of light 
nuclei (from protons to oxygen ions) in homogeneous tissue-like media. The phantom 
material and size, as well as beam parameters such as energy spread, transverse beam 
size, emittance, angular divergence, can be set via user interface commands. 

In MCHIT the energy loss and straggling of primary and secondary charged 
particles due to interaction with atomic electrons is described via a set of models 
called 'standard electromagnetic physics'. Multiple scattering due to electromagnetic 
interactions with atomic nuclei is also included in simulations. 

In each simulation step, the ionisation energy loss of a charged particle is 
calculated according to the Bethe-Bloch formula. The average excitation energy of 
the water molecule was set to 77 eV, i.e. to the value which better describes the set 
of available data on depth-dose distributions for therapeutic proton and carbon-ion 
beams. This parameter was taken 68.5 eV for PMMA, 78 eV for graphite, 86.5 eV for 
bone tissue and 70.9 eV for muscle tissue. 

Two kinds of hadronic interactions are considered in the MCHIT model: (a) 
elastic scattering of hadrons on target protons and nuclei, which dominate at low 
projectile energies, and (b) inelastic nuclear reactions induced by fast hadrons and 
nuclei (GEANT4-Documents 2006). 

The overall probability of hadronic interactions for nucleons and nuclei 
propagating through the medium depends on the total inelastic cross section for 
proton- nucleus and nucleus- nucleus collisions. Parametrised equations by Wellisch 
and Axen 1996 that best fit experimental data were used to describe the total reaction 
cross sections in nucleon-nucleus collisions. Systematics by Tripathi et al 1997 and 
Shcn et al 1989 for the total nucleus-nucleus cross sections were used for calculating 
the probability of nucleus- nucleus collisions. 

In the MCHIT model the inelastic interaction of nucleons below 20 MeV is 
simulated by means of data driven models. Above 20 MeV the exciton-based 
precompound model is invoked (Agostinelli et al 2003, Allison et al 2006). For hadrons 
and nuclei with the energies above 80 MeV, we employed the binary cascade model 
(Folger et al 2004). Exited nuclear remnants are created after the first cascade stage 
of interaction. Therefore, appropriate models for describing the de-excitation process 
have to be involved into simulation. The Weisskopf-Ewing model (Weisskopf and 
Ewing 1940) was used to describe the evaporation of nucleons from residual nuclei 
at relatively low excitation energies, below 3 MeV per nucleon. The Statistical 
Multifragmentation Model (SMM) by Bondorf et al 1995 was used at excitation 
energies above 3 MeV per nucleon to describe the multi-fragment break-up of highly- 
excited residual nuclei. The SMM includes as its part the Fermi break-up model, 
which describes the explosive decay of highly-excited light nuclei. 
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3. Time-dependent analysis of the (3 + activity 

We follow the approach of Parodi et al 2002 in calculating the time dependence of 
/3+-activity induced by therapeutic beams. As reported by Fiedler et al 2006, the time 
structure of ion beams provided by the GSI synchrotron consists of repeated particle 
extractions (spills) and pauses. It is assumed in calculations, that each beam spill has 
duration of r s with the average intensity of J (ions/s) during the beam extraction. 
The irradiation procedure consists of N spills with pauses between subsequent spills 
of t p , as given in Table [1] Both r s and r p are in the range of 1-3 s. 



Table 1. Beam parameters for 207.92 A MeV 3 He and 337.5 A MeV 12 C used 
for irradiation of graphite, water and PMMA phantoms by Fiedler et al 2006. It 
is denoted: N — number of spills , J — beam intensity during each spill, r s — 
spill duration, t p — duration of a pause between spills. 



Projectile 


Phantom 


N 


J 


Ts 






material 




(10 8 s" 1 ) 


(8) 


00 


3 He 


graphite 


120 


1.9 


1.39 


3.10 


3 He 


water 


99 


2.0 


1.35 


3.14 


3 He 


PMMA 


120 


2.0 


1.37 


3.12 


12 C 


graphite 


120 


0.9 


2.20 


2.29 


12 C 


water 


120 


0.9 


2.20 


2.29 


12 C 


PMMA 


120 


0.9 


2.19 


2.30 



The depth distributions drii(z)/dz = fi(z) of positron-emitting isotopes of species 
i along the beam axis z produced per beam particle were calculated with the MCHIT 
model. These distributions refer to the secondary nuclei at their stopping points in 
the medium. Then, the depth distribution of the i-th isotope dNi(z,t)/dz = Fi(z,t) 
during the irradiation is expressed as a function of time: 
dFi{z,t) 



Jfi(z) - XiFi(z, t) for tj-T s <t< tj , (1) 
-XiFi(z, t) for tj <t < tj + t p , (2) 



dt 

where \i is the decay constant of the i-th isotope, \ — ln(2)/T^, 2 , where Ty 2 is the 
half- life of the i-th isotope, and tj = t s + (r p + r s )(j — 1) is the time when the j-th 
spill ends, j = 1, TV. Eqs.|T|) and ^ describe the production and decay of the i-th 
isotope during the j-th spill and j-th pause, respectively. After the irradiation i-th 
isotope decays exponentially: 

Fi(z,t) =F l (z,i JV )exp(-A 4 (t-i JV )) for t > t N . (3) 
The system (HJ,© can be solved recursively (c.f Parodi et al 2002): 
Fi(z,tj) = Fi(z,tj-x) exp(-Aj(r p + t s )) 

+ F i (z,t 1 ) for j = 2, ...,7V, (4) 

Fi(z,t\) = ^1(1 - exp(-A 4 r 5 .)) . (5) 

This gives the following expression for the depth distribution of the i-th isotope at the 
end of j-th spill: 

j'-i 

F l (z,t J )=F l (z 7 t 1 )J2 e M~X l {T p + T s )n) for j = l,...,N. (6) 

n=0 
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Since the measurements of the /5 + -activity during the irradiation were performed 
only in pauses between subsequent spills, the total number of /3 + -decays per unit depth 
during the irradiation is 

rlN N ~ 1 

= £(1 - exp(-A 2 r p )) 2 FiM) ■ (7) 

i j=l 

After the irradiation the measurements were performed in the time interval from t s t 
to tfi n continuously, thus 

—J^- = ^2 Fi(z,t N )[exp(-Xi(t st -t N )) - exp(-A, (i/ in -*jv))] -(8) 

The half-life times of the isotopes included in our analysis and listed in the 
next section are much longer than the spill duration and the pause between spills: 
(t s +T p )\i <C 1. Under this condition Eqs. ([7]) and ((5|) can be simplified. The number 
of /3 + -decays per unit length during the irradiation becomes 

^1 ~ 2 jfi(z)-^[tN - y (1 - cxpi-X^))} , (9) 

while after the irradiation 

dN p+ ^ Jh{z) 

—j— ~ 2^ A . V 1 " exp(-A^Ar)) 

x [exp(-Ai(fat - ijv)) - exp(-Aj(t/i„ - t N ))] . (10) 

Here J = Jt s /(t p + r s ) is the average beam intensity calculated for the whole 
irradiation period. 



4. Comparison of numerical results with experimental data 

4-1- Depth distributions of [3 + -activity 

The spatial distributions of the /3+-activity induced by 207.92 A MeV 3 He and 337.5 A 
MeV 12 C in various phantoms were measured by Fiedler et al 2006. In order to validate 
the MCHIT model with these data we performed calculations for graphite (9x9x15 
cm 3 , p — 1.795 g cm -3 ), water with an admixture of gelatine (H66.2O33.1C0.7, 9x 9 x 30 
cm 3 , p = 1.0 g cm -3 ) and polymethyl methacrylate (PMMA, CsHgC^, 9 x 9 x 20 cm 3 , 
p = 1.18 g cm -3 ) phantoms. 

Monte Carlo calculations with the MCHIT model have provided the depth 
distributions in the phantoms for the following positron-emitting nuclei with Tin > 
10 s: 10 C (T 1/2 = 19.255 s), n C (T 1/2 = 20.39 min), 13 N (T 1/2 = 9.965 min), 14 
(T 1/2 = 1.177 min), 15 (T 1/2 = 2.04 min), 17 F (T 1/2 = 1.075 min), 18 F (T 1/2 = 109.77 
min) . It was found that these are the most abundant positron-emitting nuclei produced 
by 3 He and 12 C in graphite, water and PMMA. Much lower yields were found for 8 B, 
9 C, 12 N and 13 0. Moreover, as the latter isotopes have rather short half-life time 
T1/2 < 1 s, they decay during the beam spills and do not produce any significant 
contribution during pauses when Fiedler et al 2006 performed their measurements. 

The calculated depth distribution of /3 + -activity for 207.92 A MeV 3 He beam in 
graphite, water and PMMA are shown in Figs. [1] [2] and [3] and compared to the data 
obtained by Fiedler et al 2006 during irradiation and from 10 to 20 min after it. Since 
measured distributions were presented in arbitrary units, in Figs. [TJ [2] and [3] they 
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Figure 1. Depth distribution of /3+-activity induced by 207.92 A MeV 3 He 
beam in graphite. The distributions of /3+-decays counted during irradiation and 
from 10 to 20 min after it are shown by solid lines in top and bottom panels, 
respectively. Data by Fiedler et al 2006 are shown by points. Contributions of 
specific isotopes and depth-dose distribution are also shown, as explained on the 
legend. 



were normalised to the corresponding maxima of the calculated dNp+ / dz distributions 
during irradiation. The same weight factor was applied for plotting the experimental 
data after irradiation. 

The dNp+ j dz distributions during and after irradiation were calculated according 
to Eqs.flT]) and flSJ), respectively, The beam parameters quoted in Tabled] were used 
in calculations. The calculated dN^+fdz distributions were folded with the Gaussian 
weigth of FWHM=8 mm. This width represents the sum of (1) the average distance 
between a nucleus which emits a positron and annihilation point ~ 2 mm, (c.f. Levin 
and Hoffman 1999), and (2) a finite spatial resolution (6.5 ±2 mm) of the PET scanner 
used by Fiedler et al 2006. The shapes of the /3 + -activity distributions shown in Figs.Q] 
[2] and [3J should be compared with the corresponding depth-dose distributions (given 
in arbitrary units) in order to investigate the correlation between them. 

Figures [TJ [5] and [3] demonstrate also the dependence of the /3 + -activity 
distributions on the elemental composition of the specified target materials. In 
graphite, which contains only carbon nuclei, mostly n C and 10 C isotopes are produced 
by 3 He via the removal of one or two neutrons from target 12 C nuclei. As shown in 
Fig.[U these nuclei are evenly distributed within the range of 3 He projectiles and n C 
is the most abundant /3 + -emitter. The MCHIT model predicts a bump near the Bragg 
peak due to 14 and 13 N, which is, however, not visible in the data. A smaller bump 



PET monitoring with 3 He beams 



7 




Figure 2. Depth distribution of /3+-activity induced by 207.92 A MeV 3 He beam 
in water. The distributions of /3+-decays counted during irradiation, from 10 to 
20 min and from 120 to 130 min after it are shown by solid lines in top, middle 
and bottom panels, respectively. Data by Fiedler et al 2006 are shown by points. 
Contributions of specific isotopes and depth-dose distribution are also shown, as 
explained on the legend. 



due to 13 N is also present in the total /3 + -activity distribution calculated 10-20 min 
after irradiation. 

A larger set of isotopes is produced by 3 He in water, see Fig. [5J The most 
abundant positron-emitting nuclei are 15 and n C. While 15 is produced by the 
removal of a single neutron from a target 16 nucleus, the production mechanism of 
n C, via the 16 0( 3 He,4p4n) 11 C reaction, is more complicated. This is reflected in the 
fact that the overall shape of the total /3 + -activity distribution during irradiation is 
satisfactory reproduced by the MCHIT model since it is mostly due to 15 0. However, 
there is a big discrepancy between theory and experiment for the time interval 10-20 
min after irradiation, as shown in the middle panel of Fig. O It can be explained 
by the deficiency of the model in calculating 16 0( 3 He,4p4n) 11 C reaction which can 
proceed through various intermediate states and reaction channels. 

The MCHIT model predicts a noticeable contribution to the /3 + -activity from 14 
during irradiation of water by 3 He. As the half-life time of 1 O is much shorter than 
that of n C, 14 gives larger contribution during irradiation, while n C contribution 
dominates 10-20 min after irradiation. The model also predicts a small bump close to 
the Bragg peak due to 17 F. A similar but somewhat shifted bump is also seen in the 



PET monitoring with 3 He beams 



8 




z (mm) 

Figure 3. Same as Fig. [2] but for 207.92 A MeV 3 He beam in PMMA and for 
/3+-decays counted from 180 to 190 min after irradiation (bottom panel). 

data for the total /3 + -activity distribution. The distal slope of the activity distribution 
can be used for determination of the 3 He range in tissues, similar to the proposal by 
Parodi et al 2002 for therapeutic proton beams. 

We have found that a better way to monitor the 3 He range in tissues can be 
provided by measuring the /3 + -activity at later times, e.g. in the time window of 
120-130 min after irradiation. The activity distributions for this time interval are also 
shown in Fig. [2] Here, the long-living 18 F isotope from the 16 0( 3 He,p) 18 F reaction 
dominates, and the peak in the activity distribution clearly marks the position of the 
Bragg peak. 

As shown in Fig. [3j in PMMA two dominating /3 + -cmitters, n C and 15 0, 
are produced by 3 He. These isotopes are produced in the 12 C( 3 He,a) 11 C and 
16 0( 3 He,a) 15 reactions on carbon and oxygen nuclei from PMMA. Due to PMMA 
chemical composition, n C is more abundant than 15 in this material. The activity 
distribution in PMMA calculated for later times, e.g. 180-190 min after irradiation, 
also has a bump close to the Bragg peak. This activity peak is due to 18 F produced 
in the 16 0( 3 He,p) 18 F reaction. 

The MCHIT model is also verified with the /3 + -activity distributions measured 
by Fiedler et al 2006 for 337.5 A MeV 12 C beam in graphite, water and PMMA, as 
shown in Figs. QJ [5] and [51 In these phantom materials 10 C and 11 C can be produced by 
single or double neutron removal from both 12 C projectiles and 12 C target nuclei. As 
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a result, the /3 + -activity distribution is characterised by sharp peaks due to projectile 
fragmentation and plateau due to target fragmentation, as seen, in particular, in Fig. [5] 




°0 20 40 60 80 100 120 140 

z (mm) 

Figure 4. Depth distribution of /3 + -activity induced by 337.5 A MeV 12 C beam 
in graphite. The distributions of /3 + -decays counted during irradiation and from 
10 to 20 min after it are shown by solid lines in top and bottom panels, respectively. 
Data by Fiedler et al 2006 are shown by points. Contributions of specific isotopes 
and depth-dose distribution are also shown, as explained on the legend. 

The overall shapes of the total /3 + -activity distributions are satisfactorily 
described by the MCHIT model. However, as a rule, the peaks in theoretical 
distributions are located 5-10 mm deeper compared to the experimental ones. This 
shift is caused by the overestimation of the 11 C production in the binary cascade model 
at low energies, as explained in Section 

As seen in Figs. HI [5] and [5] for all three materials, only n C survives within 10 
min after irradiation. However, during irradiation the contributions from 14 and 15 
are important for water and PMMA. These positron-emitting oxygen isotopes are also 
produced beyond the Bragg peak by energetic secondary protons emitted in nuclear 
fragmentation of 12 C beam. 

4-2. Total yields of positron-emitting nuclei produced by 3 He and 12 C beams 

The MCHIT model can further be validated by confronting calculated total production 
yields of specific positron-emitting nuclei with the yields measured by Fiedler et al 
2006. In Tables [21 [3] and [4] we present the total yields of the most abundant isotopes, 
10 C, n C, 13 N and 15 0, produced by 3 He and 12 C beams in graphite, water and PMMA 
phantoms. The values are given in % per beam particle. 
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z (mm) 

Figure 5. Same as Fig. [4] but for by 337.5 A MeV 12 C beam in water. 

From inspecting Table [Done can conclude that the calculated yields of 11 C are 
in better agreement with experiment than 10 C yields. The yields of n C produced by 
3 He are well described, while the model overestimates the production of 10 C by both 
beams. 



Table 2. Computed number of /3+-emitters (in % per beam particle) from 
interactions of 3 He and 12 C in graphite. Data are from Fiedler et al 2006. 





3 He, 130.03A MeV 


3 He, 166.05A MeV 


3 He, 207.92A MeV 


12 C, 337.5A MeV 




MCHIT 


Experiment 


MCHIT Experiment 


MCHIT Experiment 


MCHIT 


Experiment 


10 c 


0.712 


0.44±0.07 


1.07 0.69±0.10 


1.51 0.98±0.15 


4.56 


1.88±0.28 




5.28 


6.9±1.0 


7.99 10.0±1.5 


11.09 13.9±2.1 


34.80 


24.9±3.7 


13 N 


0.508 




0.460 


0.358 


0.155 




15 o 


0.0021 




0.0022 


0.002 


0.031 





One can also compare the model with the data on 13 N and 15 production in 
irradiation of water phantoms by 3 He and 12 C. Calculated and measured yields for this 
case are presented in Table [3J The yields of 15 produced by 3 He are underpredicted 
by the model by ~ 30% for all 3 He energies, while the production of 15 and n C 
by 12 C beam is well described. On the other hand, the model completely fails in 
describing 13 N and 11 C production by 3 He, as these yields are underestimated by a 
factor of three. We attribute this problem to the complexity of 16 0( 3 He,3p3n) 13 N and 
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z (mm) 

Figure 6. Same as Fig.H but for by 337.5 A MeV 12 C beam in PMMA. 

16 0( 3 He,4p4n) 11 C reactions. The ability of the model to predict 10 C yields for 3 He 
beams depends on the beam energy. The calculations agree better with the data at 
207.92 A MeV than at 130.03 and 166.05 A MeV. 



Table 3. Computed number of /3+-emitters (in % per beam particle) from 
interactions of 3 He and 12 C in water. Data are from Fiedler et al 2006. 





3 He, 130.03A MeV 


3 He, 166.05A MeV 


3 He, 207.92A MeV 


12 C, 337.5A MeV 




MCHIT 


Experiment 


MCHIT 


Experiment 


MCHIT 


Experiment 


MCHIT 


Experiment 


10 c 


0.039 


0.21±0.03 


0.071 


0.19±0.03 


0.118 


0.18±0.04 


1.93 


0.78±0.12 


n c 


0.415 


1.85±0.28 


0.784 


2.49±0.37 


1.28 


3.23±0.48 


13.1 


12.6±1.9 


13 N 


0.154 


0.49±0.07 


0.235 


0.80±0.12 


0.321 


1.02±0.15 


1.02 


2.40±0.36 


15 o 


2.30 


4.40±0.66 


3.75 


6.29±0.94 


5.65 


8.29±1.24 


15.3 


14.6±2.2 



Calculations and experimental data for PMMA phantoms irradiated by 3 He and 
12 C are presented in Table HI The calculated yields of 10 C and 13 N produced by 3 He 
in this material are well described by the MCHIT model, while the yields of most 
abundant n C and 15 are underestimated by ~ 30%. The model is quite successful 
in describing 13 N, n C and 15 produced by 12 C beam, but the production of 10 C is 
overestimated. 

In summary, the yields of the most abundant isotopes n C and 15 produced by 
12 C in water and PMMA are very well described by the MCHIT model, see Tables |3] 
and [4j However, the yields of n C and 15 0, which are abundantly produced also by 
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Table 4. Computed number of /3+-emitters (in % per beam particle) from 
interactions of 3 He and 12 C in PMMA. Data are from Fiedler et al 2006. 





3 He, 130.03A MeV 


3 He, 166.05A MeV 


3 He, 207.92A MeV 


12 C, 337.5A MeV 




MCHIT 


Experiment 


MCHIT 


Experiment 


MCHIT 


Experiment 


MCHIT 


Experiment 


10 c 


0.398 


0.38±0.06 


0.607 


0.53±0.08 


0.862 


0.68±0.10 


3.30 


1.64±0.25 


n c 


2.89 


4.70±0.71 


4.47 


7.05±1.06 


6.26 


9.64±1.45 


23.5 


22.0±3.3 


13 N 


0.332 


0.23±0.05 


0.337 


0.28±0.06 


0.321 


0.44±0.09 


0.425 


0.63±0.13 


15 


0.87 


1.53±0.23 


1.39 


2.35±0.35 


2.06 


3.19±0.48 


5.10 


5.14±0.77 



3 He in water and PMMA, are underpredicted at all 3 He energies. This means that 
there is a room for improvement of the GEANT4 nuclear reaction models with respect 
to 3 He-induced reactions. 

5. Calculations of /3 + -activity distributions in tissues 

As demonstrated above, the total yields and spatial distributions of /3 + -activity 
produced by therapeutic beams depend essentially on the elemental compositions 
of target materials. Therefore, for studying the feasibility of the PET monitoring 
method in real tissues irradiated with 3 He we have performed calculations for two 
homogeneous phantoms with elemental composition similar to muscle (9x9x30 cm 3 , 
p = 1.061 g cm~ 3 ) and compact bone (9 x 9 x 15 cm 3 , p = 1.850 g cm~ 3 ). The 
elemental composition was taken in the following mass fractions: H - 10.2 %, C - 
14.3%, N - 3.4 %, O - 71%, Na - 0.1 %, P - 0.2%, S - 0.3 %, CI - 0.1 %, K - 0.4% for 
muscle tissue, and H - 6.4 %, C - 27.8%, N - 2.7 %, O - 41%, Mg - 0.2 %, P - 7%, S 
- 0.2 %, Ca - 14.7 % for compact bone tissue. The beam parameters in calculations 
of 3 He irradiation of muscle (bone) were taken the same as for 3 He beam in water 
(graphite), as listed in Table [TJ The calculated depth distributions of /3 + -activity in 
tissues irradiated by 207.92 A MeV 3 He are shown in Figs. [7] and [8] 

Due to the presence of Na, P, S, CI, Ca, K and Mg in muscle and bone tissues, 
19 Ne (T 1/2 = 17.22 s), 21 Na (T 1/2 = 22.49 s), and 30 P (T 1/2 = 2.498 min) can also 
be produced in fragmentation reactions in addition to the isotopes analysed in the 
previous sections. However, only negligible yields of 19 Ne and 21 Na are predicted by 
the MCHIT model and they can be safely neglected. The /3 + -activity distributions 
in muscle are similar to those in water (c.f. Fig. [2]). In bone tissue, however, 30 P 
is additionally produced by the 31 P target fragmentation. In fact, both 30 P and 15 
noticeably contribute to the dNp+ jdz after the distal edge of the Bragg peak, as shown 
in Fig. [H] This is similar to n C distribution in graphite and PMMA irradiated by 
3 He (see Figs. Q] and [3] where also a tail of the /3 + -activity is present beyond the Bragg 
peak). 

Only n C and 18 F survive in muscle and bone tissues at later time after irradiation 
The presence of 18 F opens a new way to monitor the 3 He range in tissues, as the 18 F 
peak (see the bottom panels of Figs. [7] and [5]) clearly mark the position of the Bragg 
peak. 
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Figure 7. Depth distribution of /3+-activity induced by 207.92 A MeV 3 He beam 
in muscle tissue. The distributions of /3+-decays counted during irradiation, from 
10 to 20 min and from 120 to 130 min after it are shown by solid lines in top, 
middle and bottom panels, respectively. Contributions of specific isotopes and 
depth-dose distribution are also shown, as explained on the legend. 



6. PET monitoring with low energy proton, 3 He and 12 C beams 

The nuclear pick-up reactions leading to the production of 14 0, 17 F and 18 F nuclei 
have the maximal cross sections at low energies. In this energy regime the velocities of 
projectile nuclei are comparable to the velocities of intranuclear nuclcons due to Fermi 
motion. This gives optimum conditions for transferring nucleons from one collision 
partner to another during their collision and enhances the production of 14 0, 17 F and 
18 F. 

It is instructive to consider the distributions of /3 + -activity in muscle produced by 
various low-energy beams during irradiation, as shown in Fig. [5J In these calculations 
the time structure of all three beams was assumed the same as for graphite irradiation 
by 3 He, see Table. [TJ It is expected that nuclear transfer reactions are more important 
at low energies, while nuclear fragmentation reactions contribute less because they 
have certain energy thresholds. 

As one can see in Fig. [§1 the distribution of positron emitting nuclei produced 
by low-energy 87 MeV proton beam is almost uniform and poorly correlated with the 
position of the Bragg peak. In fact, the Bragg peak in the depth-dose distribution is 
located at the region with a negligible /3 + -activity. A similar distribution is predicted 
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for 15 produced by 102 A MeV 3 He. However, 17 F nuclei are additionally produced 
by 3 He in nuclear charge pick-up reactions close to the Bragg peak. This makes the 
total distribution of /3 + -activity for 3 He more suitable for the determination of the 3 He 
range in tissues by the PET method, as the distal end of the /3 + -activity distribution 
marks clearly the position of the Bragg peak. It is advisable to perform the PET 
measurements with low-energy 3 He beams to quantify the contribution of charge pick- 
up reactions. The /3 + -activity distribution produced by 162 A MeV 12 C is also suitable 
for PET monitoring due to the presence of broad peaks associated with 10 C and 11 C 
nuclei. 

7. Reliability of calculational results 

We discuss the discrepancies between the MCHIT results and the experimental data 
by Fiedler et al 2006. The calculations agree with the data on the total yields of 10 C, 
n C, 13 N and 15 in graphite, water and PMMA at — 30 - 50% accuracy level, see 
Tables [5J31 However, the agreement with the activity distribution measured 10-20 
min after irradiation of water by 207.92 A MeV 3 He is poor, see the middle panel of 
Fig. [3j We identify this discrepancy with the poor description of 16 0( 3 He,4p4n) 11 C 
reaction, and we conclude that it has to be improved in GEANT4. In this reaction 
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Figure 9. Depth distribution of of /3+-activity during irradiation by 87 MeV 
protons (top panel), 102 A MeV 3 He (middle panel) and 162 A MeV 12 C beam 
in muscle tissue. 



a compound nucleus 19 Ne can be created leading to a larger longitudinal momentum 
transfer as compared with the direct mechanism. Therefore, the n C nuclei produced 
in decays of 19 Ne will have on average a larger longitudinal momentum and will stop 
closer to the Bragg peak. Therefore, the agreement in shapes of the calculated and 
measured distributions of 11 C nuclei can be possibly improved by taking into account 
the formation of 19 Ne. 

The model was also confronted with the measured activity distributions from 
337.5 A MeV 12 C in graphite, water and PMMA. The largest discrepancy between 
calculations and data obtained 10-20 min after irradiation was found for the graphite 
phantom, see the bottom panel of Fig.[4j To identify the origin of this discrepancy the 
cross section er( n C) of the 12 C( 12 C,n) n C reaction was calculated with the GEANT4 
toolkit. This was done following the expression: 

o-( 11 C)= d -^ n-\ (11) 

dz \ z =o 

where dNn c /dz^ z=Q is the number of n C nuclei per incident 12 C ion per unit length 
at the entrance point z = to the graphite. In fact, we have averaged dNnc/dz over 
the region < z < 1 mm. The 12 C concentration in graphite is n = 9. • 10 22 cm~ 3 . 
The obtained results are compared with the data by Yashima et al 2003, 2004 in 
Table O One can see that the MCHIT model overpredicts the cross section tr( 11 C) 
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by 30-60%. The ratio CTRg ANT4 /<7ifg seems to grow at lower beam energies. This 
explains the shift of the peak of the calculated activity distribution to larger z with 
respect to the data. 

Table 5. Cross section of the reaction 12 C( 12 C,X) n C calculated using GEANT4 
in comparison to the data from Yashima et al 2003, 2004 . The error bars on 
theoretical results are pure statistical. 



E/A, MeV 


CT GEANT4 j mb 


mb 


100 


144±14 


88.3±3.2 


230 


106±17 


79.0±7.9 


400 


100±10 


68.6±2.5 



The production of 18 F by 36-40 MeV 3 He ions in water was studied by Fitschen et 
al 1977 and by Knust and Machulla 1983 . They reported 20 mCi//xA and 19 mCi/^A 
activity of 18 F after 3 and 2.5 hours of irradiation with 3 He beam, respectively. The 
activity calculated by the MCHIT model for similar irradiation conditions amounts 
to 18.7 mCi//iA, which is in very good agreement with the experimental yields. This 
analysis shows that the low-energy nuclear data play an important role in designing 
reliable models for heavy-ion cancer therapy. 

8. Summary and conclusions 

We have considered nuclear reactions induced by 3 He and 12 C beams in tissue- like 
materials from the view point of their suitability for PET monitoring. As found, in 
addition to nuclear fragmentation reactions of projectile and target nuclei, leading to 
creation of 10 < 1:L C and 14 ' 15 0, the contributions from several nuclear pick-up reactions, 
12 C( 3 He,X) 13 N, 12 C( 3 He,n) 14 0, 16 0( 3 He,X) 17 F and 16 0( 3 He,p) 18 F must be taken 
into account. The pick-up of nucleons is quite efficient at low collision energies 
when the relative velocities of nuclei are comparable to the characteristic velocities 
of intranuclear nucleons due to their Fermi motion. On the other hand, it is known 
that this is the region of ion-beam energies with the maximum of relative biological 
effectiveness (RBE). 

As pointed out a long time ago by Osgood et al 1964, Cirilov et al 1966 
and Hahn and Ricci 1966, there exist different mechanisms of pick-up reactions. 
Indeed, the proton stripping reaction 12 C( 3 He,d) 13 N is a direct reaction, while the 
12 C( 3 He,n) 14 process goes through the formation of 15 compound nucleus. In both 
cases, the reaction cross section increases at low projectile energies. As follows from 
our calculations, 13 N and 14 are concentrated near the Bragg peak of 3 He in graphite 
and PMMA (Figs.[T]and[3]), while the distributions of 10 C, n C and 15 are flat because 
such nuclides are resulting from target fragmentation reactions. 

The process 16 0( 3 He,p) 18 F is of particular interest since 18 F is a long living 
isotope (Tj/2 = 109.77 min), making it suitable for off-line monitoring of /3 + -activity. 
According to Hahn and Ricci 1966, the cross section of this reaction has a maximum 
436 ± 44 mb at Ei a b = 6.3 MeV and drops at higher beam energies. As demonstrated 
above, 2-3 hours after irradiation by 3 He a peak in the /3 + -activity distribution due 
to 18 F is developed in water and PMMA. The position of this peak well matches the 
position of the Bragg peak. Counting statistics in the experiment by Fiedler et al 2006 
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should, in principle, allow to identify this peak. Since the measured biological wash- 
out time is about of 91-124 min in muscle (Tomitani et al 2003), i.e. comparable to 
physical half-life of 18 F, the biological wash-out should not drastically decrease PET 
signal even for water-dominated tissues. In bone tissue PET signal is expected to be 
robust and survive 2-3 hours after irradiation. Alternatively, the biological wash-out 
of /3 + -activity at later times (> 1 hour) may be studied with 18 F. 
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